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Abstract 
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for loop graphs is proved. Corrections to the Newtonian potential are also 
derived. 
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1 Introduction 



In the field-theoretical approach to general relativity, Einstein's theory can be 
reduced to a theory of gravitation in flat space-time by an expansion of the 
gravitational Lagrangian density, the Hilbert-Einstein Lagrangian density C HE , 
as an infinite series in powers of the gravitational coupling constant. 

In such an expansion the inherent non-linearity of Einstein equations ap- 
pears as a non-linear interaction between gravitons due to their gravitational 
weight. Thus, Einstein's theory will be considered as a self-interaction of or- 
dinary massless rank-2 symmetric tensor gauge fields in flat space-time. The 
Lorentz covariant quantization program for the so obtained field theory was 
proposed in Q, g, g, [§, §, [§ and @ (and references therein). 

For an ample treatment of this subject, see [||, || and p0| . 

Although covariant Feynman rules in quantum gravity (QG) were soon de- 
rived [11], [^] which allowed a calculable perturbative expansion of QG [ 13 1 , 
it was soon realized that in the standard perturbative framework radiative 
corrections within the theory were plagued by severe ultraviolet (UV) diver- 



gences [pL4 1 , fig] (for a review, see ||T6|, |L7]], (T 



It turned out that pure (that is without matter fields) QG was one-loop finite 
due to the Gauss-Bonett identity in four dimensions. The obtained one-loop 
divergences are such that they can be transformed away by a field renormal- 
ization. But two-loop calculations J19|] , | |20| | and |2lJ] yield non-renormalizable 
divergences. 

In the meantime, it was realized that the reason for the UV divergences 
lies basically in the fact that one performs mathematically ill-defined opera- 
tions, when using Feynman rules for closed loop graphs, because one multiplies 
Feynman propagators as if they were ordinary functions. 

Therefore a new strategy was developed in order to avoid the appearance 
of UV divergences once and for all. This was done by Epstein and Glaser in 
the early seventies p2| , then further applied to QED by Scharf |23| and to 
Yang-Mills theories by Diitsch et al. |24|] . 

In the resulting scheme, called 'causal perturbation theory', the central ob- 
ject is the ^-matrix, whose perturbative expansion is computed taking causality 
as cornerstone so that all expressions are finite and well-defined. UV finiteness is 
then a consequence of a deeper mathematical understanding of how loop graph 
contributions have to be calculated. 

Therefore, power-counting perturbative renormalizability no longer repre- 
sents a criterion for distinguishing viable theories from ill-defined unrealistic 
theories. 

Within the causal perturbation scheme, one-loop contributions to graviton 
self-energy are calculated in this paper and shown to be UV finite without 
the introduction of a regularization scheme and therefore cutoff-free. The ob- 
tained graviton self-energy tensor depends logarithmically on a mass scale, which 
breaks scale invariance, and satisfies the appropriate gravitational Slavnov- 
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Ward identities (T^j, p5| , |2q ], only if graviton and ghost loops are added up 
together. 

An aspect of causal theory applied to QG relies in the fact that one works 
with free graviton fields in a fixed gauge, therefore, for general gauge calculations 
we refer to (27j, f28|. 

Although in this paper we do not present two-loop calculation, the causal 
method ensures us of their UV finiteness [29]. 

For the explicit quantization of the graviton field and the subsequent con- 
struction of the physical subspace of the graviton Fock space which contains 
physical graviton states and for the proof of unitarity of the S'-matrix restricted 
to the physical subspace, we refer to |30|| , which also provides us with the basic 
notations and definitions. 

QG coupled to photon fields and to scalar matter fields within causal per- 
turbation theory is considered in Q] and [^| , respectively. 

The paper is organized as follow: in the next section, after a brief intro- 
duction to causal perturbation theory, the transition from general relativity to 
perturbative quantum gravity in the causal approach is grounded and the con- 
dition of perturbative gauge invariance is presented and some consequences are 
drawn. In Sec. || the inductive construction of the graviton self-energy is ex- 
plicitly carried out, the issues of non-normalizability of QG and distribution 
splitting are touched on. In Sec. || the Slavnov-Ward identities for the two- 
point function under investigation are verified and their relation to perturbative 
gauge invariance is discussed. In Sec. [5] the freedom in the normalization that 
is inherent to the causal inductive construction is settled. Corrections to the 
Newtonian potential due to graviton self-energy are discussed in Sec. ||, while 
perturbative gauge invariance to second order in the loop graph sector is shown 
in Sec. [?]. In the technical appendices, the formulae needed for the causal con- 
struction of the 2-point distributions and for the sum of self-energy insertions 
are derived. 

We use the unit convention: h = c = 1, Greek indices a, f3, . . . run from to 
3, whereas Latin indices i, j, . . . run from 1 to 3. 



2 Graviton Coupling and Perturbative Gauge Invari- 
ance 

2.1 S-Matrix Inductive Construction 

The central object of causal perturbation theory |23[| , p9| , |30[| is the scattering 
matrix S. Being a formal power series in the coupling constant, we consider it 
as a sum of smeared operator-valued distributions of the following form: 

°° 1 f 

S(g) = 1 + ^2— ld i xi...d i x n T n (xi,...,x n )g(x 1 )-...-g(x n ), (2.1) 
n=i n ' J 
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where g is a Schwartz test function (g E 5(R 4 )) which switches the interaction 
and provides a natural infrared cutoff in the long-range part of the interaction. 

The 5-matrix maps the asymptotically incoming free fields on the outgoing 
ones and it is possible to express the T n 's by means of free fields. In causal 
perturbation theory interacting quantum fields do not appear. 

The n-point operator-valued distribution T n is a well-defined 'renormalized' 
time-ordered product expressed in terms of Wick monomials of free fields. T n 
is constructed inductively from the first order T\(x), which corresponds to the 
interaction Lagrangian in terms of free fields, and from the lower orders Tj, 
j = 2, . . . , n — 1 by means of Poincare covariance and causality. 

Causality leads directly to UV finite and cutoff-free T n -distributions in every 
order without introducing any counterterm. 



2.2 First Order Graviton Interaction 



Following the usual approach |fL3| , [23] , we start from the Hilbert-Einstein La- 
grangian (without cosmological constant) 

C HB = ^-^R, (2.2) 

where R is the Ricci scalar and k 2 = 2>2ttG with G = Newton's constant. We 
use the same notations as in |3(| and |p3| . Expanding the Goldberg variable 
9^ u := V~9 9^ U i n an asymptotically flat geometry 

g»»(x) = V " U + Kh^{ X ), (2.3) 

where rf" v = diag(l, —1,-1,-1) is the flat space-time metric tensor, we find the 
non-terminating expansion of C HE 



oo 

C he = J2^^b, (2.4) 

j=0 



where C%^ E represents an interaction involving j + 2 gravitons. Eq. (2.3) defines 
the dynamical graviton field h^ u (x) propagating in the flat space-time geometry. 

The lowest order C%' E is quadratic in hP ,v (x) and in the Hilbert gauge 
h^ v {x)^ = the graviton field h fJ,,/ (x) obeys the wave equation 

Dh fMU {x) = 0. (2.5) 



Since the perturbative expansion for the S'-matrix (2.1) is in powers of the 
coupling constant k, we consider the normally ordered product of the first order 
term in ( ]2.4| ) 



T*(x) = i K -.C^ E (x):=i^{+ :bP a h t f p h°f: -\ :h^h p h a : + 

+ 2 : h pa h u £h pa p : + : h pa h <a h pa a : -2 : h pa h ap h a ^ : } , (2.6) 
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as the cubic interaction between gravitons or first order graviton interaction. 
For brevity, we omit the space-time dependence of the fields, if the meaning is 
clear. 

For convenience of notation, the trace of the graviton field is written as 
h = /i 7 7 and all Lorentz indices of the graviton fields are written as superscripts 
whereas the derivatives acting on the fields are written as subscripts. All indices 
occurring twice are contracted by the Minkowski metric rf" v . 

The 'non-renormalizability problem' of quantum gravity arises because of 
the presence of two derivatives on the graviton fields in ( |2.6| ) whose origin lies 
in the dimensionality of the coupling constant: [k] = mass -1 . 



2.3 Quantization of the Graviton Field, Perturbative Gauge In- 
variance and Ghost Coupling 

We consider the graviton field h fJjU ( X^j ELS cl free quantum tensor field which sat- 
isfies the wave equation ( |2.5| ) and quantize it by imposing the Lorentz covariant 
commutation rule 

[h a P(x),h» v (y)] = -ib a ^ u D (x-y), (2.7) 

with 

yxfav ._ 1 (^fc^v + _ ^/y^ . (2.8) 

and Dq(x) is the mass-zero Jordan-Pauli causal distribution: 

1 



D (x) = D { +> (x) + D { ^(x) = ^S(x 2 ) sgn(x u ) 



(2tt): 



■ Jd 4 P 5{p 2 )sgn{p )e- ip x . 



(2.9) 



The gauge content of quantum gravity is formulated by means of the Lorentz 
invariant and time-independent gauge charge Q: 

Q:= J d 3 xh a P(x) :f3 d° u a (x) , (2.10) 

x°=const 

where ^{x) is a C-number vector field satisfying Ou^(x) = 0. The gauge charge 
generates the infinitesimal operator gauge transformation of the graviton field 

dqh^ix) := [Q, h^(x)] = -i b^ U P(x), a . (2.11) 

The ten components of the symmetric rank-2 tensor h) iV contain more than the 
true physical degrees of freedom of a massless spin-2 particle, this additional 
freedom could be suppressed by a gauge condition h!% = and a trace condition 
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h = 0. As in gauge theories these conditions are disregarded at the beginning 

and considered later as conditions on the physical states. In [ 30 1 the explicit 

construction of the Fock space for the physical graviton states is carried out 

and there it is shown that the physical subspace can be defined as T p hys = 
ker({Q,Qt}). 

Gauge invariance of the S-matrix means 

lim [Q,S{g)) =0. (2.12) 

This condition is fulfilled, if the perturbative gauge invariance condition for the 
n-point distribution T n 

dgT n (xi, . . . , x n ) = sum of divergences , (2-13) 

holds, because divergences do not contribute in the adiabatic limit g — ► 1 due 
to partial integration and Gauss' theorem. 

For n = 1 the above requirement is not at all trivial, because dqT^{x) ^ 
divergence. This requires the introduction of an interaction between graviton, 
ghost and anti-ghost, i.e. the first order ghost coupling [35], [36 

!? = »«(+: u»Xp uP ■ ~ ■ u> VP <p ■ ~ ■ ^ P <9 ■ + : ~<^ V < P ■ ) • ( 2 - 14 ) 
The ghost fields must be quantized as free fermionic vector fields 

a u u (x) = 0, au u (x) = 0, {u^{x),u v {y)} = irf v D^{x - y) , (2.15) 

whereas all other anti-commutators vanish. The ghost and anti-ghost fields 
undergo the infinitesimal operator gauge variations 

d Q u a (x) := {Q,u a (x)} =0, d Q u a (x) := {Q,u a {x)} =ih^(x) 3 (2.16) 

under the action of Q, so that the sum of fl2.6| ) and ( |2.14| ) preserves perturbative 
gauge invariance ( |2.13|) to first order: 

d Q T^ u {x) = d Q {T*{x) + T?(x)) =: &* V T» /X {x) . (2.17) 



One possible form of Ty t (x), the so-called Q-vertex, was derived in |33| . 

The definition of the Q-vertex from Eq. ( [2.17 ) allows us to give a precise 
prescription on how the right side of Eq. ( |2.13| ) has to be inductively constructed. 
We define the concept of 'perturbative quantum operator gauge invariance' by 
the equation 

n d 

dgT n (xi,... ,x n ) = ^-Q-^T^^xi,... ,x h ... ,x n ), (2.18) 
i=i Xl 

where is the 'renormalized' time-ordered product, obtained according to 
the inductive causal scheme, with a Q- vertex at xi, while all other n — 1 vertices 
are ordinary T\ -vertices. 
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2.4 Consequences of Perturbative Gauge Invariance to Second 
Order for Two-Point Distributions 

We derive now some consequences from the condition of perturbative gauge 
invariance to second order for loop graphs. From the structure of 
it follows straightforwardly that by performing two field contractions ([O]) the 
resulting 2-point distribution T2{x,y) will be of the form 

T 2 (x,y) =+ :h^(x)h^(y): it hh (z) a ^+ : (x) >e h^ (y) : it dhh (z)^ u + 

+ :h a P(x)hr(y), p : i t hdh (zf^ u + : h"? (x) , e hT (y), p : it dhdh (z)^ u 
+ :v?{x)u»{y)y. it udil (z)^+ :u a {x) il3 u 5 (y): it dau {zf^ + (2.19) 
+ :u y (x) iS u tl (y) ! p: it duaa (z)]^+ :vf 1 (x)^u 0i {y)^: i t mau (z)^ , 

where z := x — y. The subscript on the numerical t-distribution denotes the 
structure of the external fields attached to them. This T 2 (x,y) describes the 
graviton self-energy and the ghost self-energy. The corresponding tensors will 
be given in Sec. 3^ and in Sec. 7.1 , respectively. 
Perturbative gauge invariance to second order 



d Q T 2 (x, y) = d^ n (x, y) + c^T^x, y) 



(2.20) 



enables us to derive a set of identities for these distributions by comparing the 
distributions attached to same external operators on both sides of ( 2.20| ). 

The left side of fl2.20|) is obtained by calculating the infinitesimal gauge 
variations of the external fields^ by means of Eqs. ( p.ll[ ) and ( 2.16| ) and we 
isolate the terms with external operators of the type :u(x)h(y) : so that 

d Q T 2 (x,y) = +b a ^ 5 \ :^{x), s h^{y):t hh {z) a ^ v + :v? {x) Ae h^ {y) : t dhh (z)f^ 



afifiv 



+ 



(2.21) 



+ :u\x)h^{y)^ T : t u9a (z)^ p rT+ :u^x\ s h^ (y), p , T : rT + ■ ■ ■ ■ 

On the other side, d^T^^ + duT^ 2 contains also operators of this type. Using 
a simplified notation which keeps track of the field type, of the derivatives and 
of the position of the z/-index which forms the divergence in ( gjgjj ), then the 
Q- vertex of ( |2.17|) reads (see for the detailed form): 



T» /X (x) :=+ :duhd*h: 


+ 


udhd'^h 


: + 


:u v dhdh: 


+ 


:duhdh v 


+ 


+ :du u hdh: 


+ 


udhdh u 


: + 


:dudhh u : 


+ 


: d^uudu : 


+ 


+ : ududu u : 


+ 


uuddu v 


: + 


:u u uddu: 








1 cIq(:uu:) = — :u{Q,u}:= — 


:uh 


and (1q(: 


uu : ) 


— : {Q, u}u 




hu : 





(2.22) 
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while the 'normal' vertex reads 

Ti :=:hdhdh: + :dudhu: + :duhdu: . (2.23) 

Then, performing two contractions between Ty^x) and Ti(y), the contributions 
in which have external operators of the type :u(x)h(y) : are 

!?/!(*, y) =+ :rf(x), 5 h^{y): t% uh {z)f v ' + :u\x)^' (y)y. t» udh {z)]\f + 
+ :u^x)h^{y): C( Z )^+ :u^x)h^ (y), p : O) 7 ^' + 
+ t uh (z)l^'+ : ^(x)^' (y), p : W*)'^ + 

+ :^(x), 5 ^'(y): t duh {z)^ v ' + :u v {x), 5 h^ \y\ p : t 9udh (z) s \^' + 
+ :u v {x) Ae h^'{y), p : t ddudh {z)^' + :u v (x) >s , e h^' (y) : W^Jf' + 
+ :^(x), 5 , e ^'( y ), p : ^(z)^'+ :uT(x) Ae ^(y): C^)^' . 

One should not forget that there exist also terms with external operators of the 
type :u{x)h(y) : coming from T^ 2 (x,y), which is inductively constructed with 
T\(x) and T",^{y). They read 

T% /2 (x,y) =+ :u'(x)h^(y): l^z)^ '+ :^{x)h^' {y)y. l u9h {z)^' + 
+ :u\x)h^\y)y. CJ*) 7 ^ + :u^x)h^(y) jP : U*) 7 |? + 
+ :^(x)^(y): C^lT' + 

+ :u-r(x)^(y)^: W^)?^ :^(x),^'(y), p : ^(^iT' 
+ :v?(x), 5 h^(y), p : l dn9h (z)f p + :v?(x) iS h^(y): l duh (z)^ ■ 

(2.25) 

Here, the numerical distributions are denoted by I. According to Eq. ( 2.2C| ), we 
have to apply d x v to T% /1 , (|2.24D , and dl to T£ /2 , ( ^25|) . After that, we gather 
the various terms according to their Lorentz structures given by the position 
of the indices and the number of derivatives acting on the external fields. We 
compare then the C-number distributions attached to the external operators: 

:u'{x) jS h^{y): 1 {x) Ae h^ {y) : , : u* (x) }S h>" (y) , p : , 
:u^x) A eh^(y), p :, :u'(x)h^(y)^ T : , : v?{x)^{y)^ T : , :u'(x)h^(y):, 
:v?(x)h^{y)y., : v?{x) A e,^ v {y), p : , :u>{x) At ^{y) : , (2.26) 
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between dqTi and d^T^ + d„Ty 2 so that we obtain the identities 

(il), ^(z)^ = | + d?n uh {z)^ v + 4(*F^ + ^ t uh (x)^ + 

+ t guh {z)^ + d% lg uh (z)J^ + 9% l 9uh (z)J^ 

v*« t ehh ( Z )^ = { + t% uh { Z )t v + ^ w*),'r + ^ *&*(*)£r + 



(13) , w^)^ = I + a? Wjjr + 4,(^) 7 |'r + a? W*),|r + 

+ ^7T ^ludh + ^ ^"Sh ( z )l|'p + ^ ^ s "' 1 

(14) , 6°^* W^J?" = { + W(T + ^ W^J" + ^ 

(15) , V UT = { + r/ pr U^F 1 "" + Ca^F^ + *T UWt 
(i7), = |+ d^tl h {z)^ + d*t uh (z)^ + d y j: h (zT'^ + dll uh {zrA , 

m, o = |+ d%tZ ah ( z y^ + d*t udh ( z )\^ + C(z)^i^ + dy p i udh (z)^ 

(i9), = |+^ t dduah (z) s ^; + 4„ 9fe (^) , 

(<io), o = |+^ t 99uh (z) s ^ + o*) jH • 



(2.27) 

These identities hold among the C-number 2-point distributions constructed in 
second order perturbation theory. Some of them have been explicitly checked 
by calculations, but there is no doubt about their validity, because in Sec. [7.2| , 
the condition ( 2.20| ) of perturbative gauge invariance to second order for loop 
graphs is proved. 
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In the case of QG coupled to photon fields (31]] and scalar matter fields | p2[ , 
these identities are less involved and from them we can derive easily the Slavnov- 
Ward identities for the 2-point connected Green function with photon and mat- 
ter loop, respectively. We will return on these identities and their relation to 
the gravitational Slavnov-Ward identities in Sec. 4. 

3 Two-Point Distribution for Graviton Self-Energy 

It is our aim in this section to apply the causal scheme to QG in order to cal- 
culate the 2-point distribution T2(x,y) which describes the graviton self-energy 
contribution. We explain step by step how T 2 (x,y) has to be constructed ac- 
cording to the general rules of the causal scheme p3[ . 

The are two important pieces in the inductive calculation that we are going 
to carry out: the first one is the calculation in momentum space of the product 
of positive/negative parts of Jordan-Pauli distributions (see App. 1 and App. 2 
for the technical details) and the second one is the causal splitting procedure 



(see Sec. |3.5|) according to the correct singular order (see Sec. p.4[) . 



3.1 Inductive Construction 



First of all, from the commutation rules (2.7) and ( |2.15| ) we compute the con- 
tractions between two field operators: 

C{h af3 {x) h^{y)} := [h aP {xf-\ h^(y) (+) ] = -i b a ^ u D { +\x - y) , 

C{u»(x) u u (y)} := K(x) ( "\ u u (y) <+) } = +irT D^(x - y) , 
C{vP{x) u v (y)} := {u^{x) ( -\u u {y) {+) } = -i^ u D^\x-y)- (3.1) 

where (±) refers to the positive/negative frequency part of the corresponding 
quantity. 

The first step in the construction of T2(x,y) consists in calculating the aux- 
iliary distributions 

R' 2 (x,y) := -Tf+«(y)Tt+«(x), A' 2 (x,y) := -T?+ u (x)T?+ u (y) (3.2) 
from these we form the causal distribution 

D 2 (x,y) := R' 2 {x,y) - A> 2 (x,y) = [l*+«(x),I**«(y)\ . (3.3) 

Causal means that the numerical part of D 2 (x,y) has support inside the light 
cone. Being T^ +U (x) a normally ordered product, we have to carry out all 
possible contractions between the two factors in Q3.2| ) using Wick's lemma. In 
this manner D 2 (x,y) contains tree contributions or scattering graphs (only one 
contraction and four external legs), loop contributions (two contractions and 
two external legs) and vacuum graph contributions (three contractions and no 
external legs). Note that, due to the presence of normal ordering, tadpole 
diagrams do not appear in causal perturbation theory. 
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3.2 Example of the Calculation 

Let us illustrate how to construct Di{x, y) by explicitly working out an example. 
We take into account only the first term in the graviton coupling T^(x) so that, 
from the A 2 (x, y)-distribution, that can be decomposed as a sum of 25 different 
contributions Yli j=i A%(x, 2/) Cm \ we pick up only the term A' 2 {x, and 
in addition we carry out the loop generating double-contractions only between 
graviton fields that carry derivatives, so that for A' 2 (x, y) (M) we obtain 



A' 2 (x,y)^ 



— ) :h a Ph p °h p %::h^h^W s u : 



2 contractions 



+ ^ :h^(x)h^{y): 



+ C{h^(x), a h^(y)^} ■ C{h^(x),p W s {y) tV ) + 



+ C 



{h^(x), a h? 5 (y), v ) ■ C{h^(x). fi W 5 (y) 4l ] 



+ other contractions . 



(3.4) 



Using the relations in (|3.1|), we find 



A' 2 (x,yY^ = +^- :h^(x)h^(y): 



x 



y)) 



x ( _ i \r^d x p diD { +\x - y )) + {-i r^flsw 



x(-ib^ s dpyD^(x-y)) 



+ other contractions . 



Since bP^ s b pay8 = W and 3^D^\x - y) = -^D^\x - y) we obtain 



A' 2 (x,y)^ =:h<#{x)hr(y): a' 2 (x - y)%^ + 



where 



5 K 



D {+ \ n (x 



D «l\8> ~ V) ■= ^D£\x - y) ■ 8ffiD$>{x - y) . 



p,x n(+) 



y) + D ( :J l(¥ (x - y) 



y)) 



(3.5) 



(3.6) 



(3.7) 



The products between derivatives of Jordan-Pauli distributions are calculated 
in App. 1. Analogously, by taking into account that Dq(jj — x) = — Dq (x— y), 
we find 



R' 2 (x,y)^ = :h^(x)h^(y): r' 2 (x - y)™ + 



where 



r' 2 (* ~ Vtpl ■= ~ (D^ u (x -y) + D^(x - y)) 



D *t\aM ~ V) ■= d«d*D l ^(x - y) ■ dffiD^(x - y) 



(3.8) 



(3.9) 
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Therefore, according to Eq. fl3.3|), the D2(x, y)-distribution has the form 

D 2 (x,y)™ =:h^(x)h^(y): d 2 (x - y)^ + . . . , 

d2(x - ytg v = r' 2 {x - y)^ - a' 2 (x - y)%^ . (3.10) 

The most important property of D 2 (x, y) is causality, but only the numerical 
distribution d 2 (x — y) is responsible for this property: 



supp(c/ 2 (z)) C V+(z)\JV-(z), with z:=x-y, (3.11) 

(see below). The products of Jordan-Pauli distributions appearing in ( 3.10[ ) are 
easily expressed in momentum space, see App. 2, so that we obtain 

«2(F)S, = -f(ri3 F e(p 2 )e(V), 

^(p)S& = ©b 2 ) ©(V) ; (3.12) 

and therefore 

*(P)S& = P (P)%, U ®(P 2 ) [©(P°) - ©(V)] = ^ F ©(P 2 ) sgn(p°) , 

(3.13) 

where Pip)^^ is a Lorentz covariant polynomial of degree four, this degree is 
given by the number of derivatives on the contracted lines. Causality is evident 
from the scalar distribution d(p) := @(p 2 )sgn(p°). For z 2 < 0, we may choose a 
reference frame in which z a = (0, z), so that 

d(z) = jdp S gn(p°) jd 3 P e(p 2 -p 2 )e +i P- z = 0, (3.14) 

because of the signum- function in p°. Therefore d(z) vanishes outside the light 



cone, see Eq. (3.11). 



3.3 Causal D 2 (x, y) -Distribution for Graviton Self-Energy 

The total D 2 (x, y)-distribution for the graviton self-energy through a graviton 
loop is obtained by calculating the 25 contributions coming from the graviton 
coupling Ti, not only the terms with two external graviton fields without deriva- 
tives, but also these with one or two derivatives. In addition, there are also 16 
contributions coming from the ghost-graviton coupling where one performs two 
ghost-anti-ghost contractions. Summing graviton loop and ghost-anti-ghost 
loop contributions we obtain ^ 



D 2 (x,y) =+ :h^{x)h^(y): df {x - y)'^ 



+ 



+ :h^(x), 7 h^(y): d™ (x - y)^ v + 
+ :h^(x)h^(y), p : d* b \x-y)%^ + 
+ :h^(x) >7 h^(y\ p : i$>{x - y)^ . 



(3.15) 



2 the notation '•[•' keeps track of the exact position of the indices 
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The tensorial distributions have in momentum space the structure 

d?(p)ip flv = P ii, (p)lp fi J(p), i = 4,3a,36,2. (3.16) 

The results for the tensorial distributions, being too long, are not given here. 
See App. 4 for the explicit form of the distributions appearing in Eq. ( 3.15| ). 



These will be used to calculate the graviton self-energy tensor in Sec. |3.6 . 

In order to obtain T2(x,y), we have to split the /^-distribution into a re- 
tarded part, i?2j and an advanced part, A2, with respect to the coincidence 
point z := x — y = 0, so that supp(i?2(^)) Q V + (z) and supp(yl2(.z)) Q V~(z). 
The correct treatment of this coincidence point constitutes the key to control 
the UV behaviour of the 2-point distribution. 

This splitting procedure affects only the numerical distributions d 2 in J3.15Q 
and must be accomplished according to the correct singular order u;(4 l> ) of the 
distribution. It describes the behaviour of d 2 (z) near the coincidence point 
z = 0, or that of d 2 (p) for p — ► 00. If ui(d 2 ^) < 0, then the splitting is trivial. 
On the other side, if uj(d 2 ) > 0, then the splitting is non-trivial and non-unique: 



4 l '(^ - V)tp,u — 4\x - y)^ + £ {Ca, t D*}^ 8V(z - y) , (3.17) 

|a|=0 

and a retarded part r 2 (x — y) is best obtained in momentum space by means of 
a dispersion-like integral, see Sec. |3.5| , which, however, presents some difficulties 
in the massless case. 

The C a /s in Eq. (jTlTD are undetermined but finite normalization constants 



and D a is a partial differential operator acting on the local distribution 5 (4) (x — 
y). The second term on the right side of Eq. ( [3.17 ) represents therefore a freedom 



in the normalization which is inherent to the causal splitting of distributions and 
will be discussed in Sec. || by taking physical conditions into account. 

In the case of Eq. ( [3.16 ), we find from direct inspection of the distributions 



that 

w(4 4) )=4, w(4 3a) ) = 3, w(4 3i,) ) = 3, w(4 2) ) = 2. (3.18) 

The singular order depends on the structure of the graph, namely on the number 
of derivatives acting on the contracted internal lines of the loop. For a precise 
formulation, see below. 

The last step in the inductive construction of T 2 (x, y) consists in subtracting 
R' 2 (x, y) from B,2(x, y), see Sec. |3.5| . The singular order remains unchanged after 
distribution splitting: cj(c?2) = uj{r 2 ) = a; ^2). 

3.4 Singular Order in Quantum Gravity 



Before undertaking the splitting of D 2 (x, y) according to Eq. ( 3.17 ), we give the 
formula for the singular order of arbitrary n-point distributions in perturbative 
quantum gravity. 
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We consider in the n-th order of perturbation theory an arbitrary n-point 
distribution T£(xi, . . . , x n ), appearing in Eq. (|2.1|) , as a sum of normally ordered 
products of free field operators multiplied by numerical distributions 



n h n u n a 

T°(xi,... ,x n ) =: Y[h(x kj ) J\u(x m ) JJn(a? ni ): t°(xi,... ,x n ). (3.19) 

j=l i=l 1=1 

This T° corresponds to a graph G with external graviton lines , n u external 
ghost lines and riu external anti- ghost lines. The singular order of G then reads 

uj(G) < 4 — nh — n u — riu — d + n . (3.20) 

Here d is the number of derivatives on the external field operators in ( 3. 19| ) . 
The '<' means that in certain cases the singular order is lowered by peculiar 
conditions, e.g. by the equations of motions of the free fields. 

The explicit presence of the order of perturbation theory renders the theory 
'non-normalizable', that is the theory has a weaker predictive power but it is 
still well-defined in the sense of UV finiteness. 

We give some hints of the inductive proof of ( 3.20|) |23|| , |37f . First of all, 
the assumption ( |3.20| ) must be verified for n = 1: uj(T^ +u (x)) has to be zero 
(because of ui(5(x)^ = 0), a result which is correctly given by ( 3.20| ) after direct 
inspection. 

In the inductive construction of T n from the T m 7 s, m < n — 1, we must 
consider tensor products of two distributions 

T rj i(xi,... ,x r )T s>2 (yi,... ,y s ), (3.21) 

with known singular order w(T ri i) = uj\ < 4 — n/ ll — n ui — n„ 1 — d\ + r and 
^(T s ^) = 0J2 < 4 — rih 2 — n U2 — nn 2 — d 2 + s. According to the inductive 
construction, this product has to be normally ordered giving origin to all possible 
contraction configurations. We assume that / contractions arise during this 
process. Taking translation invariance into account the numerical distribution 
of the contracted expression is of the form 

I 

ti(xx -x r ,... ,av_i -x r ) d a J D { +) (x r . -y Sj )t 2 (yi ~y s ,--- ,Vs-i ~ Vs) = 
= t(£i,... ,^ r _i,r/i,... , 77s— i , rj), (3.22) 

with £j := Xj—x r , rjj := yj — y s , rj := x r —y s , aj = 0, 1, 2 and a = Y?j=i a j- Then, 
using the distributional definition of the singular order [23], we may conclude 
that 

u (t) =UJl + bj2 + 2l-A + a. (3.23) 
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Inserting the expressions for w% and uj 2 in Eq. (3.23), we get 



lo{1) < 4 - (n hl + n h , 2 + n Ul + n U2 + n Sl + n Sa - 21) - (d\ + d 2 - a) + (r + s) . 

(3.24) 

The first bracket represents the number + n u + ria of external fields after I 
contractions, the second bracket gives the number d of derivatives remaining on 
these external fields, if the I contractions carry a derivatives. Since r + s = n, 
Eq. (ET20D is proved. 

In the usual QFT formulation, Eq. ( 3.2C| ) would imply that QG is 'non- 



renormalizable', since uj(G) increases without bound for higher orders in the 
perturbative expansion. This means that there is a 'proliferation' of divergences 
and of counterterms to compensate them. 

The situation is different in causal perturbation theory: we are facing in this 
case a 'non-normalizable' theory, i.e. each of its diagrams is finite due to the 
causal splitting method, but the number of the free, undetermined and finite 
normalization constants in ( |3,17D increases with n. The question is then to find 
enough physical conditions or requirements to fix this increasing ambiguity in 
the normalization. 

3.5 Splitting of the D 2 (x, ^-Distribution 

We now carry out the splitting of the distribution Di{x,y) in Eq. ( 3.15| ). 



Let us consider for example the numerical tensorial distribution d^ip)^ 



which has singular order four from Eq. ( 3.1 §| ) of from Eq. ( |3.20| ). Because 



of the decomposition ( 3.16j ), it suffices to split the scalar distribution d{jp) 



0(p 2 )sgn(p°) with io(d) = and then multiply the so obtained retarded part by 

the same tensor -P (4) (p)i/3^^ as given by E q. (|3.16| ). 

Usually, a special retarded part in Eq. ( |3.17D , if it exists, is given in momen- 
tum space by the dispersion integral 



which is called 'central splitting solution', because the subtraction point |g3[ 
is the origin. But this formula cannot be used directly in the case of massless 
theories, because the integral is divergent. In order to circumvent this deficiency, 
we shift the original distribution d(x): 

d q (x) := e iq - x d(x) , d q {p) = d{p + q), q 2 < , (3.26) 



so that the central splitting solution r~ (p) of the shifted distribution exists 

We cannot obtain the retarded part of the original distribution simply by 
letting q — > 0, so we take advantage of the local ambiguity in the splitting 
procedure and consider another retarded part of d g (p), given by r q {p). Since 
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two retarded distributions differ only by local terms in configuration space, we 
obtain in momentum space for fixed q that the difference reads 



r q {p)-r°Jp) = P q (p), 



(3.27) 



where P q (p) is a q-dependent polynomial in p of degree lo. Then we construct 
r(p), a retarded part of the original distribution d(p), from ( 3.27] ) by taking the 
limit 



r(p) = lim f q {jp) = lim [rg(p) + P q (p)] ■ 



(3.28) 



Here the addition of the g-dependent polynomial P q {p) must be accomplished in 
such a way that the limit exists. This corresponds to a finite renormalization. 

Using ( 3.25 ) with ( |3.26| ) in ( |3.28| ), we obtain for p G V + , with q — ► in such 
a way that p — q G V + , q 2 < 0: 



r(p) = lim 



lim 



i f +0O i 

— / dt- 

-oo 



d q (tp) 



lim 



— / dt 

2vr 



2^ 



t - tO) (l - t + iO) 
d(tp + g) 
_ x v t - iO) (1 - t + iO) 



(3.29) 



(* - iO) (1 - 1 + iO) 



e((tp + g) 2 ) sgn(tp° + ?°) + P,(p) 



The zeros of (tp + g) 2 are ti^ 



p ■ q ± VN) with N := (p ■ q) 2 - p 2 q 2 , so 



that the integral in (|3,29|) may be simplified to 



*2<0 

^ I dt + — 



2tt 



2tt 



di J- ( - 



ti>0 



1 



t- 1 



t7rtf(t-i; 



2^ 



i 7T + log 



+ O0 



and the limit reads 



r(p) = lim 

g^O 



27 



? 7T + log 



+ o(v1?|)) +p,(p) 



(3.30) 



(3.31) 



Being u = 0, we can add the polynomial P q {p) := 27 log {\Q 2 \/M 2 ) , where M > 
is an arbitrary mass scale, so that we obtain a Lorentz invariant splitting solution 



f(p) 



i{ log {W 2 



I 7T 



P eV + . 



By analytic continuation in M 4 + iV + , we find for pGl 4 



f( P r = — 10, 



-p 2 - ip u 
M2 



(3.32) 



(3.33) 
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As pointed out at the end of Sec. |3,3] , the T2(x, y)-distribution is obtained 
from R2{x,y) by subtracting R' 2 (x,y). This subtraction affects only the scalar 
distributions. Since f'(p) = —0(p 2 )<d(—p°), we obtain 

t(p) = f(p) an - r'(p) 

= i( log (^) W)©(p 2 )) +©(p 2 )e(-p°) {3U) 
= i f log ( tig) - i = i ior f - p2 - z0 



2vr V b \ M 2 J ^ ') 2vr fa V M * 

The ambiguity in the normalization present in the splitting D 2 — > R2 + N 2 , 
Eq. ( |3.17| ), will be discussed in Sec. || 

3.6 Graviton Self-Energy Tensor from the T 2 (a?, ^-Distribution 

Gathering all the results of the previous sections, Eqs. ( |3.15| ), ( |3.16| ) and ( 3.34| ), 
we find that the 2-point distribution that contributes to the graviton self-energy 
reads 



T 2 (x,y) =+ :h^(x)h^(y): t^(x - y)^ u + 



+ :h^(x)^(y): t f a \x - y^„ + 



(3.35) 

+ :h^(x)h^(y), p : tf\x-y)%^ + 
+ :h^(x)^(y), p : 1$>(x - y)^ , 

where the tensorial distributions have in momentum space the structure 

if(p)'Uu = p{l) iP)'Uu Kp) , i = 4, 3a, 36, 2 , (3.36) 

where i(p) is given in Eq. (3.34) and the polynomials are those of Eq. (3.16). 

The ^-distributions appearing in Eq. ( 3.35| ) have already been introduced 
in Eq. fl2.19| ): t { 2 4) = it hh , t { 2 a) = it ahh , t ( 2 b) = it hah and t { 2 ] = it ah8h . 

Since divergences in the adiabatic limit of Eq. ( |2.1[ ) do not contribute, we can 
obtain from T2(x,y) by partial integration the graviton self-energy contribution 

T 2 (x,y) hSE =:h^(x)h^(y): iU(x - y\^ v . (3.37) 

The main result of our calculation in second order causal perturbation theory 
is the graviton self-energy tensor IL(x — y) a p P u which is given by the following 
combination of t2(x — ^-distributions 

iU(x - yU, v := + tf{x - y)+ - d^\x - y)%, v + 

1 1 (0.00) 

-1-0 z 2 {x y) af3fll/ o o i 2 {x y) a/3fiu , 
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where we have carried the derivatives acting on the external fields in Eq. ( |3.35| ) 
on the corresponding ti{x — y)-distributions. In momentum space, it reads 



n(pL 



n(p)^ loop +nb)f;;: loop 

- 656 p a p f3 p fl p u - 



k 2 it 



960(2tt) e 



208 p 2 {p a p li rf v + p»p u r, af) ) + 



+ 162 p 2 {p a p^if v + pVj]* + + ^P v r, a ^) + 



162 p 4 (ri a ^ u + ri au V ^) 



+ 118;pW 1/ 



log 



-p 2 - iO 
M 2 



(3.39) 



Separate calculations for the graviton loop and ghost loop give the following 
contributions to the graviton self-energy tensor, respectively (see App. 4): 



grav. loop 



fr^Nghost loop 



880, -260, +160, -170, +110 



log ((-p 2 - iO)/M 2 ) , 



+ 224, +52, +2, 



log((-p 2 -iO)/M 2 ) ; (3.40) 



where we have adopted the convention of writing only the coefficients of the 
tensor according to the structure given in Eq. (|3.39|) and H := K 2 7r/960(27r) 5 . 
Our result, numerical coefficients and logarithmic dependence on p 2 /M 2 , 

obtained using ad- 



agrees with the finite part of previous calculations [25|, 
hoc regularization schemes. As a consequence, the absence of UV divergences 
means that we do not need to add counterterms [O], [15| involving four deriva- 
tives to the original Hilbert-Einstein Lagrangian in order to obtain UV finite 
radiative corrections to the graviton propagator. In our approach, all the ex- 
pressions are cutoff-free and finite at each stage of the calculation due to the 
causal scheme. 



4 Gravitational Slavnov— Ward Identities and Pertur- 
bative Gauge Invariance 

Gravitational Slavnov-Ward Identities 

The gravitational Slavnov-Ward identities (SWI) p5| , [26|, [38|, [39| are derived 
in standard quantum field theory from the connected Green functions. We 
construct the 2-point connected Green function as 

G(p)Z,u ■= W Dob) n^) 7 ^ W Do(P) > 

where DK[p) = (2ir)~ 2 (— p 2 — iO) . The two attached lines represent free 
graviton Feynman propagators. Then the gravitational SWI reads 

P a P"G( P ) l X v = 0. (4.2) 
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Since the tensorial structure of a general self-energy tensor may be character- 
ized by the five coefficients A, B, C, E, F in the standard representation tl{p) = 
E [A, B, C, E, F] log ((-p 2 - iO)/M 2 ) as in Eq. ( ggg ), then the SWI are equiv- 
alent to the following relations 

A 

— - B + E + F = 0, C + E = 0. (4.3) 
These relations are satisfied by the coefficients of the self-energy tensor in 



Eq. ( 3.39 ), only if both graviton and ghost loops are taken into account. There- 



fore our result satisfies the SWI. 

SWI and Other First Order Couplings 
If we had chosen another ghost coupling instead of the one in Eq. (2.14), for 



example T?(x) = : u^h^u p p : +2 : u p ph v Pu% : ) @, then we would 

have violated the SWI, because in this case the new graviton self-energy tensor 
through graviton and ghost loop would have had the form 

fl(p)^ = H[-836,-238,+^,-162,+118] log (( -p 2 - iO)/M 2 ) , (4.4) 

so that it would have not satisfied the SWI in Eq. (|4.3j ). 

Analogously, if we had disregarded the last two terms in the graviton cou- 
pling T±(x), Eq. ( pH|) , being a divergence due to the presence of two equal 
derivatives, then we would have violated the SWI, too, because in this case we 
would have obtained a 'reduced' graviton self-energy tensor through graviton 
loop of the form 

ftO)^ = 3 [ - 880, +160, -108, +58, -62] log ((-p 2 - iO)/M 2 ) , (4.5) 

so that the sum of graviton and ghost loop would have not satisfied the SWI in 
Eq. ([4,3j). This happens although the difference between the two tensors can be 
written as a divergence 

n(p)^ lo ° P - n(p) r ^, = d*W(x -y) = divergence , (4.6) 

due to the vanishing of the coefficient proportional to p OL pl 3 p^p v and therefore 
this difference should not be physically relevant in the adiabatic limit g — > 1 of 

S(g)- 

SWI and Perturbative Gauge Invariance 

The SWI for QG coupled to photon fields ]3l| and matter fields |32| are equiv- 
alent to the transversality of G(p)^ , namely p a G(p)^ . In the case of 
self-coupled and ghost-coupled gravitons, the complexity of the gauge structure 



implies that only the weaker condition (4.2) can be satisfied by the graviton 
self-energy tensor, which is not transverse. 
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Perturbative gauge invariance to second order for loop graphs can be formu- 
lated by means of the identities (il), ... , (ilO) in Eq. Q2.27|) . Note that these 
identities can be easily generalized to the n-th order for graphs with two exter- 
nal operators ('legs'). The identities will imply a relation among the graviton 
self-energy tensor 

u( P r^ = t hh (py^ + i Pe t Shh (p)^ - i Pp t hah {p) a ^ v + PePp t dhah (p)f p ^ , 

(4.7) 

the ghost self-energy tensors t um (p) and t 8uda (p) and distributions coming from 
the Q-vertex, namely t uh , t u9h and others with one graviton and one ghost 
as external legs. We work in momentum space, but we skip the hat on the 
distributions. 

From the structure of the first four identities in ( 2.27| ), it follows straight- 
forwardly that 



irp s p T t u8h (pyr + 7 



+ 



+ P 2 ldu8h(p)J\ + 



(4.8) 



+ PePX iLahipVsir + VePv lau8h(p)i + 7 <-> & + M ^ v 



Here we have used (il), ■ ■ ■ , (i4) to express the graviton self-energy distribu- 
tions by means of the Q- vertex distributions and also (i9), (ilO) to simplify 
the expression. Symmetrization of the left side is indeed necessary, because the 
right side is symmetric. With (id), we can write the second bracket as 



PePr 



+ rf T t au9a (p)]^ + n <-> v + 7 



(4.9) 



The connected 2-point Green function has a second 6-tensor attached on the 
right, therefore we compute 

+ \ V lS tu h {pt a + \ ^ t„ h (p) kp - \ rf 5 U{p)^ rf° + 7 <- 5 + 

- ^4,>) 7 i'r - |pe4,(p) 7 i'r + j^wtv CT +7~ s+ 



+ 



(rT + M <-» v + 7 - 5)2 



(4.10) 
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We apply now j> 7 to the above expression. This enables us to use other identities 
Ui7) and (*8)) , so that after a little work we get 



i PSPe tudh (p) + "~ PSPe t u8 h (P) V 



pp pa 



+ 



PePrP-, 



+ 



+ 



A ( m \S\po 



- Pe luenip) 5 ^ ~ Pe La h (p) SWp ~ Px & h (p) 
~P\ lu dh (p) 5 \° P ~ P° ^dhipft ~ Pp iuahipft 



W 7 {PeLM^+PxltMT + PxL Sh (p)°\ 



A 



,S\X 



(4.11) 



Using now the identity (i5), we can express the /„ eh -distributions by means of 
the t^-distributions: 

p 7 (b^ n{p) a ^ v b^ p °) = l - + Pl t 5 uh ( P y\p° - ip 7 t s uh { P y iW v pa 



+ Pst uh {p)\ ptJ - l -P5t uh (pt^ prj 



-iPSPe t u9h (p)\ P e a + T^fWe iu8h(p)jg' i V 

+ P<j iauec(p)jj e P + P P *a«ss(p)j|f - Pi/ £a«e«(p)jj* ^ + 7 ^ # 
- P<t tuaaCp)*'^ - tuduip) \e + Pv tnduip) 5 ^ V pa 



p,p pa 



PeP-y 



iPe 
2 



+ 



(4.12) 



This expression is the desired identity that involves the graviton self-energy 
tensor n(p) Q/3/ ^, the ghost self-energy tensors t 9u9a ( P ) and t u9a ( P ) and the dis- 
tributions coming from the Q-vertex. Schematically: 



p 7 (bU(p) by° pa = X(t uh (p),t u9h (p), . . . )*>!» + (t 9um (p) + t uX ( P )Y' pa , (4.13) 

where X S,P<J represents the first eight distributions in ( [4.12 ). In QG, (bH( P ) b) is 
not transversal, because the right side of ( 4.131) does not vanish. Note that the 
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remaining Q-vertex distributions cannot be eliminated from (4.13) as already 
noticed in pq ], although in a very different approach to gauge invariance in QG. 
By multiplying ( 4,12; ) with p p , we obtain 



PpP~i 



7|W 



1 



+ P8Ppt uh {p) lp(7 ~ ^P6Pat uh { P yW 



+ 



PeP-yP 



2 r 



tdudu{p)j\ e ^~t du8il (p)\ 



i PpP-fPe t u9h (p) |'; + ^ PvPlPe Kah (P) 

i PSPe t udh (p) jf + - PSPe t udh (p) 
o2 



7|W 



ip e p 



tudu{p) 



5\<J 



(4.14) 



This equation does not a priori imply the SWI (4.2). We can summarize 
Eq. flOD as 

p p Pj (bll(p) by 5pa = p p X(t uh (p),t udh (p), . . . ) 5,pa +p p (t du9a (p) + t uea (p)) 5 ' pa . 

(4.15) 

With the explicit results for t dudtl (p) and t uail (p) obtained in Sec. 7.1, we find 
that there must be a compensation between Q-vertex contributions p p X s,pa 



and ghost self-energy contributions so that the right side of ( 4.15 ) vanishes 
identically, because this has already been verified by explicitly calculating the 



left side of (4.14). This may also explain why in the two cases with different T" 
or different T±, investigated above, the SWI of Eq. (4.2) or Eq. (fD^) are not 



fulfilled: Eq. ( |4.15| ) still holds true, therefore the theory is still gauge invariant 
to second order, but both sides do not vanish identically. 



5 Fixing of the Freedom in the Normalization of the 
Self-Energy Tensor 



As seen in Sec. 3.4, in causal perturbation theory the problem of eliminating 
infinitely many UV divergent expressions is changed into the issue of fixing an 
increasing number of free undetermined local normalization terms that arise as 
a consequence of distribution splitting in each order of perturbation theory. 



5.1 Normalization Terms of T 2 (x,y) 



For simplicity, we consider only the freedom in the normalization of Eq. ( 3.37j ) 
instead of Eq. ( |3.35| ). Since the singular order is four, local normalization terms 
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N%(x,y) ~ 5 (4) (x — y) of singular order 0, . . . ,4 appear during the process of 
distribution splitting and in momentum space they can be written as 



N 2 (x,y) hSE =:h^(x)h^(y): iN(d x , d y ) a ^ v 5^(x - y) , 



(5.1) 



where the odd terms are excluded owing to parity. N{py^p is a polynomial in 
p of degree % with % = 0, 2, 4. We assume in addition that only scalar constants 
should be considered, because vector-valued or tensor-valued constants may 
enter in conflict with Lorentz covariance. Therefore we make the following 
ansatz taking also the symmetries of Tl(x — y) a /3/j,v into account 



N(p) (2) 



N(p) (4) 



+ C 4 (p a PiMV/3u + PaPvV0n + PpPfiVav + PfiPv^) + (5.2) 
+ C 5 P 2 (Va^Pv + VavVPv) + C6 p 2 ^a/?^ , 
C7,C 8 ,C 9 ,Cio,Cn 



ci, . . . ,cn are undetermined real numbers. Requiring the SWI to hold, we can 
reduce the normalization polynomials to 



N(p) m 



N(p) (2) 



N(p) (4) 



(2) 

a/3fj,u 
(4) 

a/3/j.u 



0. 



0,c 5 + c 6 ,-c 5 ,c 5 ,c 6 l p 

C7, Cio + Cn + — , -Cio, Cio, en 



(5.3) 



in such a way that only five undetermined coefficients remain to be fixed. The 
self-energy tensor supplemented by the normalization terms then reads 



n(p)^ = n(p)^ + JVfr)^ + N( P )^ U . 



(5.4) 



5.2 Total Graviton Propagator and Mass and Coupling Con- 
stant Normalization 

The task of eliminating the remaining freedom in Eq. ( |5.4| ) can be accomplished 
in our case by considering the total graviton propagator as a sum of the free 
graviton Feynman propagator p0| 



(n\T{h^(x)h^(y)}\n) = -ib a ^D F (x - y) , 



(5.5) 
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from Eq. (2.7) and the contraction (3.1), with an increasing number of self- 
energy insertions. In momentum space we therefore obtain the series 



D(p)toT 



+ b a ^ 5 £>g(p) n(p)» Sp „ Dq (p) n( P )f TAK b^ b F {p) + 
DM r 



+b a ^u( P r P ^ 5 D( P ) 



I tot 



(5.6) 



where II {pY c 



b a /3tT, we find 



(27r) 4 if-(p)a8uv After multiplying with (Z)q(p)) 1 and with 



7<5 



i erpa n(p) 



pcr-yS 



ETflV ■ 



(5.7) 



Since b af 3 pu b^ u pa = l a p pa := {r) ap r}p v + rjauVp^) / 2 represents the unity for rank- 



4 tensors, from Eq. (5.7) we can derive the form of the inverse of the total 
propagator: 

(5.8) 



= (27T) 2 [b aPflu (-p 2 - iO) - (27T) 2 n(p)^ 

Our aim is to impose mass and coupling constant normalization conditions on 
the total graviton propagator, therefore we have to invert the expression in (|5.§|). 
For this purpose we express the inverse of the total graviton propagator in the 
'projection operator' Qipy^^v Das i s > introduced in App. 3 and find 



8=1 

x { + x\ (p 2 - iO) - (2tt) 2 5 xj log ((V - iO)/M 2 ) } , (5.9) 



with 



i=l 



n(p) 



i=i 



so that, with Eqs. (C.10), (C.ll), we obtain 



D{p) 



tot 
afifiu 



(2tt) 2 L 2 



+ 7} [llocfiVPv + VauVfo) 



+ 



2 Va^u _ p 2_ iQ + Mp2) 



(5.10) 



-p 2 -iO + fi(j> 2 ) 
+ non-contributing terms . (5.11) 
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The non-contributing terms vanish between conserved energy-momentum mat- 
ter tensors or between physical graviton states |30t| . 

Mass normalization (the graviton mass remains zero after the radiative cor- 
rection due to the self-energy) requires 



Mp 2 



p 2 =0 



0. 



1,2. 



(5.12) 



and coupling constant normalization (the coupling constant is not changed by 
the radiative correction) requires 



p- 



0. 



J = 1,2. 



(5.13) 



p 2 =0 



We work out these two conditions for the case j = 1. Since fi(p 2 ) has the form 
(see Eqs. Q, Q, ( pTlCD , da8|) and ( |CTlD ) 



/i(p 2 ) = -2(2vr) 2 ~ (-162p 4 l O g((V-i0)/M 2 ) +c 5 p 2 + Cl0 p 4 j , (5.14) 

the two conditions (|5.12|) and ( 5.13 ) are satisfied, if we set C5 = 0. The analysis 
for the case j = 2 is much more involved due to the complexity of Eqs. (C.8) 
and (p. 11 ), although not conceptually difficult, and it yields the condition cq = 
0. The remaining normalization constants are not determined by Eqs. ( 5.12j ) 
and ( 5.13; ), so that we are left with 

n(p)^ = H[ - 656, -208, +162, -162, +118] log + ^ 4) (pU^, 

(5.15) 



where N w (jp) a p^ v is a general SWI-invariant polynomial given in Eq. (5.3) with 
three free normalization constants. In order to simplify the above expression 
and to reduce further the freedom in the normalization, we rescale the constants 
07,010 and en in the following way 

Cll = +1181og(M 2 /K 2 ), c 7 = -6561og(M 2 /L 2 ), c w = -162 log(M 2 /iV 2 ) ; 

(5.16) 



so that, with (H ) 



2\104 



(N 2 ) S1 (L ; 2 f 2 /{K 2 )° 9 , we obtain 



n(p)S 



656 log 



~P 



iO 



L 2 



-208 log 



-p 2 - iO 



+ 162 log 



-p 2 — iO 
~Jq2 



,-162 log 



-p 2 — iO 
~N~ 2 



H 2 
118 log 



iO 



i^ 2 



(5.17) 



with the four arbitrary positive masses L, H, N and K. Now, since the splitting 
of the mass zero distribution d(jp) = 0(p 2 )sgn(j>°), requires the introduction of 
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a scale invariance breaking mass, it is natural to assume this mass scale to be 
unique, say Mq, which may correspond in case to the Planck mass, so that the 
graviton self-energy tensor including its normalization now reads 



n(p) 



JV 

a/3fj,u 



656, -208, +162, -162, +118 



log 



-p 2 — iO 



(5.18) 



Therefore, the whole ambiguity in the normalization of T2(x,y) hSE may be re- 
duced to the single parameter Mq which remains present in the theory. 



6 Corrections to the Newtonian Potential 



In the last years a group of papers appeared |41[], J42|] and |43j, based on the 
proposal of treating perturbative quantum gravity as a low energy effective quan- 
tum theory. In this approach, leading quantum corrections to the Newtonian 
potential were reliably calculated in the long range, low energy limit. 

In this section we will show that causal perturbation theory for quantum 
gravity yields the same result for the leading corrections to the Newtonian po- 
tential for heavy spinless particles described by the scalar field ip. 

The Newtonian potential between two masses m\,m,2 



V(r) 



G 



(6.1) 



can be obtained in the static non-relativistic limit of a single graviton exchange 

(6.2) 



tree diagram [41], [44| to lowest order in G, calculated from 



T?(x)=i :C M \x):=i^ :h^(x)b a(3flu T^(x):, 
where the conserved energy-momentum matter tensor reads 
T^(x) = i P (xy^(xy»-r ] ^CV(x), 



(6.3) 



and the expansion in powers of the coupling constant of the matter Lagrangian 
density has the form 



_. oo 



) 

M ■ 



(6.4) 



From Eq. (|6.2| ), we carry out the inductive construction of the 2-point distri- 
bution which describes the (pip — > ip(p scattering through one-graviton exchange 
and find 



ib a ^D F {x-y) 



(6.5) 
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The static non-relativistic limit of the right side of (|6.5|) yields (6.1) using Tm 



(2ir)5^ 5^°m and the fact that the Fourier transform of the momentum transfer 
p 2 reads (47rr) _1 . For a deeper understanding of the connection between S- 
matrix and potential, see |}I5|| . Inserting one self-energy contribution ( |5.18| ) on 
the graviton contraction of Eq. (|6.5|) , we obtain the order k 4 correction 



T 2 (x,yy™ = i^ :I%(x)I%fy): 0(x - »)^ , (6.6) 
where, with the self-energy tensor of ( |5,18 ), we have 



n(p)_? 

p 4 



n( P )^ a = 2 ^q^ . (6.7) 



Only the two last terms in the self-energy tensor will contribute in the following, 
because the others vanish when paired with conserved matter energy-momentum 
tensors. Therefore, even if we had not assumed a unique mass scale (as done at 
the end of Sec |5.2j ), the parameters L, H and N in ( |5.17| ) would not be physically 
observable. From ( |6.5D and (|6.7j ) we obtain in the static non-relativistic limit 
the genuine quantum correction to the Newtonian potential (inserting back the 
physical constants H and c) 

raim 2 ( ^ , Gh 206 1 
c 3 7r 30 r 



Vr (r) = _ ^( 1+ ™— ]. (6.8) 



The central piece in the calculation is the distributional Fourier transform of 
log (p 2 /Mg ) which yields (— 2-7rr 3 ) -1 and the Mo-dependence disappears from 
the non-local part of the final result, being proportional to <5 (3) (:e). The relevant 
length scale appearing in ( |6.8[ ) is the Planck length l p = y/Gh/c 3 . 

Our result agrees with the corresponding one in |4^| , although this represents 
only a partial correction to the Newtonian potential, because we have taken 
into account only the graviton self-energy contribution and not the complete 
set of diagrams of order k 4 contributing to these corrections, as, for example, 
the vertex correction or the double scattering. Therefore we cannot make any 



statement on the absolute sign of the correction in Eq. (|6.8|) . See also [46| for 
the corrections to the Newtonian potential in the framework of i? 2 -theories. 
A remark about different choices of or the graviton self-energy tensor 



of Eq. (i.4) leads to the same corrections, whereas that of Eq. ([4.5|), taking also 



the ghost loop contribution into account, leads to different corrections. 



7 Perturbative Gauge Invariance to Second Order 
for Loop Contributions 

For the sake of completeness, before proving perturbative gauge invariance to 
second order in the loop graph sector, we calculate the 2-point ghost self-energy 
contributions. 
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7.1 Ghost Self-Energy 

We follow the inductive causal construction described in Sec. |3| in order to 
construct the ghost self-energy contribution in second order perturbation theory. 
Starting with the ghost coupling T"(x) given in Eq. ( 2.14| ), we perform one 
graviton and one ghost-anti-ghost contraction, Eq. (|3.1|) in order to obtain the 
corresponding D2(x,y) distribution. The distribution splitting solution is the 
same as in Sec. 3J3, namely Eq. ( |3.34| ) with the singular order given by Eq. ( 3.20| ) 
in Sec. 3.4. Therefore, after all these steps, we obtain 



T 2 (x,y) =+:u^(x)u a (y)^: (.r 
+ :u a (x) t pu y (y) : t^" i,r 



+ :ui(x), /3 u a (y), u : t 



U,(3a) 

2 

M,(36) 
U,(2a) 



v). 







+ 



u 2 

+ :u a (x)^(y)^: t^ b \x-y) a . 



7« 

& v 
7 a 

v (3 



(x - y)j ~ + 



(7.1) 



The numerical distributions appearing here have already been introduced in 
Eq. (2.19) with the notations: t% = it udn , t^ 3 ^ 



it;, 



,u,(26) 

h 



i u, (2a) 

ii ^2 



it, 



dudi 



and 



tl 



it m9u . The results for these ghost self-energy distributions are: 

' (3a) (p) 7 |a/3 = *2' (3i,) (p)a/3|7 = E [~ ^PaPpP-y + 60pVy?7a/3 ~ 20p 2 PaVp-y + 

- 20p 2 p /3 7 ]aj ] log ((-/ - iO)/M 2 ) , (7.2) 



and 



t 



u,(2a) 



(P) 



7/3 1 a.v 



-t 



u,(26) 



(P) 



: .H [ + 160p 2 r) y pr) au + 240 p a p 7 <r) u p + 



+ 240 p u p 7 r]a(3 ~ 2A0p aPu r] 7 p - 4Q0 Pj ppr) au ] log ((-p 2 - iO)/M 2 ) . (7.3) 

By disregarding divergence terms after partial integration, we can recast Eq. 
(7.1) into the form 

T 2 (x,y) 9SE =:u 1 {x)u a (y): ill a (x - y)-y a + :u a (x)u'{y) : ill b (x - y) aj 

with the ghost self-energy tensors 

iU a (x - y), a := + dffi™{x - y) 7 / - d$d*t u 2 ^(x - y)f J , 
iU b {x - y) ai : 
that read in momentum space 



dp u / 3b \x-y) 



P 

a 7 



dpW 2b \x-y) 



v P. 

a 7 i 



(7.4) 



(7.5) 



n„(p) 



70 



A0p 2 p a p 1 + 20p 4 r? Q7 



log 



-p 2 — iO 
M 2 



-n*(p) 



(7.6) 



We do not discuss here in detail the normalization of this 2-point distribution. 
But an investigation analogous to that of Sec. 5.2 for the sum of the series 
with an increasing number of ghost self-energy insertions let us assume that the 
normalization terms with singular order smaller than four are set equal to zero, 
whereas those with singular order four may be absorbed in the scale invariance 
breaking mass M. 
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7.2 Perturbative Gauge Invariance to Second Order for Loop 
Graphs 



The introduction of the Q-vertex T^ix), Eq. ( 2.17] ), enables us to formulate 



perturbative gauge invariance by means of Eq. ( 2.18| ). We call a theory gauge 
invariant to second order, if T2(x,y) satisfies 

d Q T 2 (x,y) = &F£ /1 (x,y) + d»T£ /2 (x,y)^ (7.7) 

where T^ 1 (a;,y) (T^ 2 (x,y)) is the two-point 'renormalized' time-ordered pro- 
duct obtained by means of the inductive causal scheme with a Q-vertex at x (y) 
and a normal vertex at y (x). 

Since R' 2 (x,y) is trivially gauge invariant due to ( |3.2| ) and ( 2.171 ), ^ suf- 



fices to prove ( [r.7| ) with the retarded parts R2, R2/1 anc ^ ^2/2' i ns tead of the 
corresponding T-distributions. 

Taking Eq. Q and ( p7|) into account, we find that the /^-distribution 



is trivially invariant 

d Q D 2 (x,y) = d%D» n {x,y) + dlD» /2 {x,y) , (7.8) 



with the definitions 



D^(x,y):= [T^x),^)] , 

D» /2 {x,y):= [T x {x),T^{y)\ . UM) 



The question is whether an equation similar to (7.8) holds for the retarded 



parts R2, Ryi an d ^2/2 °^ ^2/1 anc ^ ^2/2' respectively. In the inductive 
causal construction, the splitting of distributions in D2, D 2 ^ and D^ 2 may give 
rise to local normalization terms, if the singular order is positive. We consider 
here only loop graphs in second order. Therefore, we must show that 

d Q R 2 (x, y) lo °P s + d Q N 2 (x, y) lo °P s = + d^R^x, y) loops + d^R^x, y) loops + 

+ d^N^x, y) lo °P s + d»N$ /2 (x, y) lo °P s 

(7.10) 

can be satisfied by a suitable choice of the free constants in the normalization 
terms N^ ^ , Ny° oph and N 2 j 2 ° ps of the general splitting solution of D 2 oops , 

JJ2/1 &n( * ^2/2 ' respectively. 

Since every splitting solution agrees with the original distribution on the 
forward light cone F+\{0}, gauge invariance ( 7.1C| ) can be violated by local 



terms with support x = y only. Hence, the crucial point is the correct treatment 
of the local terms appearing in ( [7.10 ). 



A careful analysis shows that that the only local terms appearing in ( 7,10| ) 
are those belonging to the normalizations N^ ^, N 2 ^° ops and N 2 j 2 ° ps of the 
distribution splitting. 
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Let us analyze the different terms in ( |7.10| ). First of all, the normalization 
terms N 2 (x, y) loo P s were already discussed in Sec. and in Sec |7, 1| , below 
Eq. ( |7.6[ ). They can be consistently set equal to zero. 

The gauge variation of R 2 (x, y) loo P s generates no local terms: the field oper- 
ators in the normal products undergo the infinitesimal gauge variations (2.11) 
and ( |2.16| ) and the numerical distributions remain unchanged. Taking an exam- 
ple from Eq. ( 3.35| ) 

d Q R 2 (x,y) lo ^ s = d Q [ :h^(x)h^(y): r$\x - y) a ^ u + . . .) 

=:uP{x)^(y): [ - i r$>(x - y) a ^] + . . . . 

r^ix — y)a/3iiv does not contain local terms, the same holds for b a ^ pa {x — 

Now we investigate the R^^x, y) loops -terms. Due to the great calculation 
complexity, we work out only a representative example, that still contains the 
main features. Let us choose in Ty^x) ~ uhh + uuu (from Q) the term — f : 

u'(x) a h(x)h(x)^ : and in T? +V (y), Q, the term -if : h°* \y)h{y) >a h(y) ,p :. 
We construct D^^x, y) lo °P s by carrying out two graviton contractions 

D^x, y)^ = i g (+ :u\x) n h^{y) : ( - 32 [ + Z3«, - DjJJ (s - y)) + 

+ :^(x) j7 My),a : ( + 8[ + D'J - D« + 0« - D<g] (x - y))) + . . . , 

(7.12) 



(7.13) 



which can be written in the form 

DZ^yf*** =+ :u\x) n h^{y): d^x - y) a(3 + 

+ :u 1 (x) n h(y) ya : d^x - y) a + . . . , 

with the Z) (± |' (p)-functions given in App. 1: 

#/! (PU = -2i « 2 [ + ^ - DjJ,] (p) , 

*/i(P) a = +4 [ + ^S* " ^S* + " ^ ^ • (7 - i4) 

Using the results of App. 1 and App. 2, we obtain 



d2/i(p)a0 = 24(2^ [ 2paP ^ ~ P 2 P* r l(? + P 2 P/3Va + pVl«(3] @(p 2 )sgn(p°), 
4x(pr = l>V] ©(P 2 ) sgn(p°) . (7.15) 
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In order to obtain R 2 J° ops , we split the distributions in ( 7.15 ). With ( |3.33j ) and 

assuming that the splitting of the massless distribution d(p) generates only one 
mass scale M, we find 

f$lM<>f> = 2 4(2ir)5 [ the salne 88 in log I TP 

f ^r = ^[pV]'o g (^g^), (7.16) 



so that 

R^ 2/l (x,yf ^ s =+ :u^x) n h^{y): r^(x - y) a p + 
+ :u 1 (x) n h(y) ya : r^x - y) a + . . . 



(7.17) 



The local normalization terms of the splitting are included in Ny^x, y) loops and 
reads: 

iVy^y) 100 ^ =+ :v?{x) n h^{y): n^(x - y) a/3 + 

+ :u 7 (x), 7 %), Q : n% fl {x - y) a + . . . . 



Here, the ri2/i-distributions contain free normalization constants. The terms 
R% /2 (x,y) lo °v s and NgJx, y) loops coi 
can be analogously calculated with x 
now to R 2 ^(x, y) loo P s ) we obtain 



Ry 2 (x, y) loo P s and N^, 2 (x, y) loo P s coming from the second commutator in ([/ 
can be analogously calculated with x <-> y and an overall sign change. Applying 



+ :u y (x) nfl h al3 (y): r^(x - y) a p+ :u 7 (x) >7At /i(y) iQ : r% /x (x - y) a + 

+ :u 7 (x), 7 /i^(y): r^(x - y) a p+ :^(x), 7 %), a : ^< A (x-y) a . (7.19) 

The first two addends do not contain local terms, because the derivatives act 
only on the external fields and the numerical distributions do not contain local 
terms by construction. The derivative acting on the numerical distributions in 
the last two terms changes them into 

-» P» r£ ;i (p) a p = 24(2 ^ )5 L + 2p PaPp + P r] af3 \ log ^ — 2 1 , 

-^/ifr) = 8^ bV] log (=^=^) • (7-20) 

Since these terms are proportional to the logarithms, none of them is local. The 
local anomaly producing mechanism described in [B3] does not apply in the case 
of loop graphs. This is a consequence of the presence in loop graphs of prod- 
ucts of Jordan-Pauli distributions that yield non-local terms after distribution 
splitting. 
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We turn now to the local terms in d^N^ 1 (x, y) lo °P s . Since oj(d 2 jf) = 3 and 
w (^2/i) = ^' these terms have the general form 

h^ipy* = +a lP > M V al3 + a 2P a rj^ + a&Prf" + o^pV + 

+ P 2 (a 5 ^r/ a/3 + a 6 p V + arpV") , (7.21) 



2/ 



= +a 8 r?^ + aflpV + aiopV 



where ai, . . . , aio are unknown parameters. Multiplying (|7.21j) with leads to 
P» ^/i(p) a/S = +ai P V /? + (a 2 + asJpV + («4 + a 6 + a 7 )p W + a5 pV , 

Pfi «2/l(P) = + a 8y a + ( a 9 + aio)]5 2 p a • 

(7.22) 



Therefore, the only local terms in ( [7.10 ) are those coming from ( 7.22j ). In our 



simplified example, perturbative gauge invariance then requires 

oi = a 2 + 03 = «4 + clq + aj = a$ = as = ag + a w = , (7.23) 

Obviously, these conditions are fulfilled by choosing a\ = . . . = aio = 0. One 
may convince oneself that the above example can be generalized to all second 
order loop graph contributions, because these follow the same pattern, although 
much more involved conditions as (|7.23| ) would appear. The important point is 
that d*R% /:L (x,y) loops + 0^R£ /2 (x,y) loops does not generate local terms. This is 
in contrast to the tree graph sector, investigated in [33]| . 

Let us add a final remark: if we do not fix the normalization N 2 (x, y) loo P s 
of R 2 (x,y) loops (as done in Sec. p[), then the condition of perturbative gauge 
invariance to second order forces us to choose the normalization constants c, of 
N 2 {x,y) loops and aj of Ng /1+2 (x, y) loops in such a way that 

d Q N 2 (x, y) loops = +d*N% /1 (x, y) loops + d»N£ /2 (x, y) loops (7.24) 

holds. Since ( |7.10| ) always holds among non-local terms by construction, a trivial 
solution Cj = dj = Vi, j always exists. This concludes the proof of perturbative 
gauge invariance to second order for loop graphs. 

Acknowledgements 

I would like to thank Prof. G. Scharf for his continuous and patient support, 
Adrian Miiller and Mark Wellmann for stimulating discussions. 

Appendix 1: The Z)^|^(p)-Functions 

The product of two Dq -distributions is well-defined in momentum space, be- 
cause the intersection of the supports of the two D^\p) is a compact set. The 
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product D\'(x) := Dq \x) ■ D (x) goes over into a convolution of the Fourier 
transforms Dq (p) using Eq. ([2.9[): 



1 



(2vr)2 

-1 
(2^)8 

-1 

Wf 



<fxD ± \x)-D^\x)e ip - x 

d A pi d A P2 5{p\) 8(±p°) 5(p|) 9(±p°) • ^d 4 * e 

d 4 /t<5((p - A;) 2 ) 9( ± (p° - A; )) 5(fc 2 ) 9(±£; ) 



-ix-(pi+p 2 -p) 



If derivatives are acting on the I2q (x)-distributions, then we define 



(A.l) 



(A.2) 



and so on if a different combination of derivatives acts on the distributions. 
Following the same calculation as in (|A.1|) , we obtain in momentum space: 

d%{p) = ^ Jd*k6{( P - kf) e(± (p° - k )) 5(k 2 ) e(±k°) 

x Pakp - k a k p , 



d 4 k5((p - k) 2 ) 6( ± (p° - k )) 5{k 2 ) @(±k° 



P\^ F > (2tt) 4 

x ~\~ PaPpk^ku Pakf}k^k v Ppk a k^k u -\- k a k^k^k v . (A. 3) 
Therefore we see that we have to deal with integrals of the type 

I {±) (p)-/a/afS/a^/a^ ■= J^k 5 ((p - A;) 2 ) G ( ± (p° - k°)) 5{k 2 ) &(±k°) 

x [l, k a , k a kp, kakpkft, k a kfjk^kiij , 

(A.4) 



which are calculated in App. 2. For the two examples in Eq. (A. 3), we find the 
relations 

-1 



dT(p) 



+i r 



D {± 1 (p) 



-1 



(2tt)4 L 



+ PaI {±) {p)p-I {± \pU\, 
+ PaPfi I {±) (P)^u ~ Pa I {±) (P)pnu 
-Ppl {±) {p)a^ + I {±) {p)a^u] , 



(A.5) 



between the D , (p)-functions and the I {±) {p)... -integrals. 
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Appendix 2: The J (±) (p) .. -Integrals 



For the case of Eq. ( ]A.4| ), the momenta k and p are restricted to the 

space-time regions {k 2 = 0, k > 0} and {(p — k) 2 = 0,j>o — > 0}, due 
to the 8- and ©-distributions in the integrand. Then p — k and k are time- 
like and therefore p is time-like. We choose a Lorentz reference frame with 
p a = (p ,0), p° > 0, then 



(+) 



(p ) = Jd 4 k5(p 2 - 2 Po k ) G( Po - ko) @{k ) 



5(k -\k\)+5(k + \k\) 
2 E k 

(B.l) 



with Ek = ko = \k\, so that 

IM( P0 ) = 4vr [ +O °dk [ + °°d\k\ \k\ 2 8(2 P0 (^ - k )) Q( P0 - k ) S(k( ' N) 

J-oo JO 1 



2|fc| 

" + OO 



dk° / |fc| 5(— — fc ) e(^ — A;o) 5(fc — 

-oo Jo 2 



VoJo 2 / 2 

(B.2) 

can be calculated analogously and the result in an arbitrary Lorentz 
reference frame is 

/(±)( p ) = ^e(p 2 )e(V). (b.3) 

Computing I l - ±) (p) a for p a = (po>0), po > 0, we have a non-vanishing con- 
tribution only for a = 0. We obtain for I (±) (po)o an additional factor ko in the 
integrand of flB.2| ), which is set equal to |fc|, because of the distribution 5(ko — \k\) 
and finally is set equal to po/2, because of the distribution 5(\k\ — ^). This 
leads to 

i (±) (p) a = lp a e(p 2 )e(±p ), (b.4) 

in an arbitrary Lorentz reference frame. 

For I^(p) a 8, we consider the covariant decomposition 

I {±) ( P U = A^(p 2 )p a p p + B^(p 2 ) Val3 . (B.5) 

It follows from J (±) (p) a a = (because of the factor k 2 5{k 2 ) in the integrand of 
Eq. (Q), that 

£ (± V) = V) • (B.6) 
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Then 

/ (±) (pWpV = |4 (±) (pV- (B.7) 

Calculating I (+) (p) a /3P a P^ for p Q = (po,0), po > Oj through the integral defini- 
tion ( |A.4p , an additional factor (po^o) 2 appears in the integrand, therefore we 
obtain in an arbitrary Lorentz frame 

I (±) (pUp a p P = \p A 6(p 2 ) 6(±p°) ■ (B.8) 

Comparing (jB^) with (gj) we find ^ (±) (p 2 ) = f 9(p 2 ) @(±p°), so that 

/ (±) (pW = I (paPp ~ j ^) 6(p 2 ) 9(±p°) . (B.9) 

For I (±) (p) a /3 M , if we calculate I (+) (p) a /3fiP a P l3 P fl for p a = (po,0), p > 0, we 
get a factor (po^o) 3 m the integrand, so that in an arbitrary Lorentz frame we 
have 

i (±) (pWVV = ^p 6 &(p 2 ) e(±p°) • (b.io) 

The covariant decomposition of i" (±) (p)a/3/i reads 

/ (±) CpW = C (± V)pW + D (±) (p 2 ) (paVfr+PpVcv+PrVafi) ■ (B.ll) 

Since I (±) {p) a / = 0, we obtain D (±) {p 2 ) = ^C (±) (p 2 ). On the other side, 
contracting the covariant decomposition of I (±) (p) a /3 ll with p a p^p^, we find 

/ (±) (pWpVV = y C^V) • (B-12) 
Comparing ( pl^ ) with ( |BlC| ) we conclude that C (±) (p 2 ) = f G(p 2 ) 9(±p°) and 

6 



/(±) (P)* = | ( + PaP/3P M - y (p<* ^ + Pf3 Van + Pfi Vafi) j @{±P°) ■ 

(B.13) 



We repeat this calculation scheme also for I (±) (p) a pfW) which has the covari- 
ant decomposition 



I (±) (p)af3tMV = E (±) (p 2 ) PaP(3P»P V +F (±) (p 2 ) [+PaPp Vlxv+PaP^Vpv+PaPv VPn + 
+ PpPn Vau + PpPu Van + PtiPv Vap) + G {±) (P 2 ) ( + Va\xVPv + VauV/3u + VapVn") • 

(B.14) 
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From J (±) (p) a a „j, = and 1^ = 0, we obtain 



Computing I i+) (p) a /3fiuP a P^p' J 'P u , for p a = (po,0), Po > 0, we get a factor 
(Po^o) 4 i n the integrand, so that in an arbitrary Lorentz frame we have 

i (±) (p)a^uP a P^ P u = ^ P S e( P 2 ) e(±p°) . (B.i6) 

On the other side, contracting the covariant decomposition of l^(p) a 0uv with 
papPpHpV anc j us i n g ( p3 ■ 1 5| ) , we obtain 

/ (±) (pWpVW' = j^P 8 E^(p 2 ) . (B.17) 

Comparing (pl7| ) with ( gig ), we find £ (±) (p 2 ) = ^e(p 2 )6(±p°), that implies 

with (pip F(±V) = ^e(p 2 )6(±p°) and G (±) (p 2 ) = ^9(p 2 )G(±p°) so 
that 

7T / p 2 ( 

I {±) {p)aP^u = Jo I + PaP P P ^Pv ~ ~g { + %i/ + ^ + 

+ PaPu VPn + PpPfi Vau + PfiPv Tjan + PpPv Vaf3^J + 
+ 48 ( + 7]a ' i ' npv + 7]a " 7] ^ + VapV/j^f) @(P 2 ) 0(±P°) • (B.18) 



Appendix 3: The Q (i) (p) a ^ 5/XI ,-Projection Operators 

The aim of this Appendix is to find a representation basis for rank-4 tensors, 
which allows to compute the inverse of the total graviton propagator ( |5.8| ) in 
Sec. || Let us define as in [47 1 

7 k^ky kpk v 

with 

d^iydyp , dpyCyp , &pV&Vp ^ pp ' (^'^) 

Then the so-called 'projection' operators are defined by 

Q( k )ap,nv = \ (da^e^ + d av ef3p + dppeay + dp v e a p) , 

1 2 

Q(^)a/3,pu = -j[d a pdpy + d av dfip - -d a/3 dp U ) , 
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and the so-called 'transfer' operators are defined by 



with the relations 



Q( k ) a8.au Q( k ) U) 



fiv,pa 



'Q( k )%, P a 



if i = 1,2,3,4; 
if i = 5,6. 



(C.4) 



(C.5) 



These and other relations can be easily calculated usin g (|C.2| ). We consider a 
rank-4 tensor in the standard basis as in Eq. ( 1Q9D or (M0|) by giving its five 
coefficients and disregarding the logarithmic dependence on k 2 /M 2 : 

[A,B,C,E,F](k) aPuu . (C.6) 

We rescale it by dividing it by k 4 , so that we can express the obtained 'rescaled' 
tensor T(k) a p^ v in the projector basis given by ( |C.3| ) and (C.4): 



T(k) 



3=1 



U) 
B,p,v • 



(C.7) 



Comparing (|C.q ) with ( |C.7| ), we find the relations between the xj coefficients 
and the A, . . . ,F coefficients: 

Xl = 2(C + E), x 2 = 2E, x 3 = 2E + 3F, 
x 4 = A + 2B + 4C + 2E + F, x 5 = x 6 = V3(B + F) . (C.8) 

The inverse of T{k) a p lil , in Eq. ( |C.7|) satisfies 

(n k )) a ^(nk)-r"=i a ^, (c.9) 



being l a g pcF = {rj afJL r]p v + rja^rjpfj) /2 the unity for rank-4 tensors and it is given 

by 



(TO*)- 1 ) 



adfiu 



+ fQ (4) -|Q (5) -fQ (6) )(fcW 



(CIO) 



=:X>Q(fc)< 

3=1 

where A := x 3 ■ X4 — £5 ■ xq. The proof of Eq. ( C.lOj ) simply consists in carrying 
out the product in ( C.9 ) and using the relations of Eq. ( C.5 ). Returning back 
to the standard representation and multiplying with k , we obtain 

2y 2 



(T(fc)- 1 ) 



a/3fj,u 



2?/i + 



2/3 



2y 5 y 5 y 3 y 2 



yi _ m m m _ m_ 
2 2 ' 2 ' 3 3 



(C.ll) 



{ k )a3f±i/ 
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Using the definitions of the y,'s as a functions of the Xj's from (C.10) and the 
inverses of the relations in ( C.8 ), we can then find the coefficients of the inverse 
(T(/c) _1 ) q ^ i/ in terms of the original coefficients A, . . . ,F. 

Appendix 4: Numerical Distributions in Eq. ( |3.15|) 

In this appendix we present explicitly the numerical distributions appearing in 
Eq. ( 3.15Q or in Eq. ( [3. 16 ). We separate the various contributions according to 
their singular order and to the presence of graviton or ghost loops. 

For the graviton loop, the distribution with singular order four reads 



df(p) 



k 2 it 



960(2vr) 4 



80, +60, +10, -70,0 



Q(p 2 ) S gn(p°) 



(D.l) 



The distribution with singular order two can be written as 



df{p) 



i\p 



K 2 TT 



960(2vr) 4 



io A (p) 



a/3fiu 



0(p 2 ) sgn(p°) , 

(D.2) 



with 



(a 



0) + (p v) 
(D.3) 



+ { PaPll (16^7?^ + 12 77/77/ + lOfy'V) + ( 

+ (a <-> f3, p <-> f)} 
+ pV { - 34 77^77^ + 60 % M ry^„ + 60 r/^r/^ } 
+ {4p Q p 7 ( - 77^77/ - r}p u r}P + 77/v) + (a «-> /3) 

+ (a <-> p, 7 <-> p, [5 <-> v) + (a <-> i/, 7 <-> p, /3 <-» p)} 
+ {p«P p ( - 3077^77^ - 3077^77^ + 1877^77^) + (a /3) 

+ (a <-> /x, 7 <-> p, /3 <-> 1/) + (a <-> i/, 7 <-> p, /? <-» p)} , 



and 



p ii(p)c$lv = + Vaf3V/j.vV 7p - 877^(77^77^ + (p <-> z^)} 

- 3 {r7a,3 W + (M <- ")} - 3 {wJvf + (D.4) 
+ { 2 rtaVfaVJ 3 + 5 Vap-VSv/ + (p ^ v) + (a ^ (3) + (a ^ /3, p, ^ v)} . 

This distribution yields after distribution splitting and after having multiplied it 
with p 7 Pp the following contribution to the graviton self-energy through graviton 
loop: 



P 1 p P if{p)^ v 



K 2 TT 



960(2vr) 4 



1200, -380, +450, -520, +330 



t{p) . (D.5) 
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The distribution with singular order three attached to :h al3 (x)^h fiu (y): reads 



°2 VPjapfiv 



iK 2 lT 

960(2vr) 4 



with 



and 



QrWoSiu, = ~ 120 V&v {Pa V + (« <- /?)} 

+ 20 {VafSPfiPuP 1 + V^PaPpPv + riJPaPpPv) , 



(D.7) 



Qn(p)apixu =+P 1 { + H0^^a/3 ~ 210 77^7?^ - 210 1/^1/^ } 

+ { Pa (-20 + 130 + 13077^,) + (a <- /?)} (D.8) 

+ {p M (-20 + 20^7?/ + 207/Jr/^) + (/x i/)} . 

This distribution yields after distribution splitting and after having multiplied it 
with i the following contribution to the graviton self-energy through graviton 
loop: 



iPltt^P)* 1 ^- 960 ( 2vr )4 



7l 



K 2 7T 



+ 200, +30, -150, +210, -110 



t{p) . (D.9) 



The distribution d^ip) a ^ v with singular order three attached to the operator 
: h a P \x)h^ v '(y) iP : is the same to the previous one with a global sign change 
and the replacements p, a, v <-> (3 and p <-> 7. After distribution splitting 
and multiplication with — ip p , it gives the same contribution to the graviton 
self-energy through graviton loop as the previous one. 
For the ghost loop we have: with singular order four 



K 2 7T 



+ 64, -8, +12, +8, +8 



e(p 2 )sgn(p°), (D.10) 



^ v 960(2vr) 4 

after symmetrization in (a/3) <-> (pv); with singular order three we have 



IK 2 TT 



+ 20p 2 {paVwVf + ia" P)} 
+10 p 2 { Pfl (riaprjf - 7] u(3 r] a p - 77^777/) + (p <-> v)} G(p 2 ) sgn(p°) , 



(D.ll) 



whereas the distribution {p)ap\j,v w ith singular order three attached to the 
operator : h a/3 (x) a hy- v (y) : is the same to the previous one with a global sign 
change and the replacements p <-> a, v <-> (3 and p <-> 7. With singular order 
two we obtain 



^ 2) iv> 7|p 

"2 yPiafiiiv 



K 2 TT 



+ 7// t/^ (80p aPll + 40 p 2 r?^)] 9(p 2 ) sgn(p°) , 

(D.12) 



960(2tt) 4 

which must be symmetrized in (a <-> /3) and (/x <-> v) and then one applies p^p p - 
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